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Abstract. The observational data on the large scale structure (LSS) of the Universe 
are used to establish the upper limit on the neutrino content marginalized over all 
other cosmological parameters within the class of adiabatic inflationary models. It is 
shown that the upper 2a limit on the neutrino content can be expressed in the form 
Vt^h'^/N^-^^ < 0.042 or, via the neutrino mass, < 4.0eV. 



INTRODUCTION 

In the light of recent data on the matter content in the Universe (0.25 < Qm < 
0.5) coming from cosmological observations and lower limit on the neutrino con- 
tribution {Qi, > 0.001) which follows from the Super-Kamiokande experiment, the 
question whether cosmology still allows neutrino mass of a few eV is of interest 
for Cosmology as well as Particle Physics. Relict massive neutrinos as collisionless 
particles with substantial velocities suppress the power spectrum of matter density 
perturbations below of free-streaming scale, which depends mainly on the neutrino 
mass, the Hubble constant and time. Hence, the comparison of the growth of mass 
fluctuations on these scales can be translated into an upper limit for the neutrino 
mass. From the Ly-a forest in quasar spectra in together with a conservative imple- 
mentation of other cosmological constraints a limit of rrii, < 5.5eV at 95% CL. has 
been reported [6] . The problem of this determination is that the damping of small 
scale power due to neutrinos can be imitated by variations of other cosmological 
parameters like the spectral index Ug, the cosmological constant A, spatial curva- 
ture Qk, the Hubble parameter h (= ifo/(100km/s/Mpc)) and the baryon content 
Qb- Therefore, Qi, (or m^) must be determined simultaneously with all these pa- 
rameters using a wide range scale of cosmological observations. Such investigations 
show also that mixed dark matter models with cosmological constant (AMDM) can 
explain virtually all cosmological measurements [14]. The goal of this paper is to 
determine more carefully the upper limit on the neutrino content and its rest mass 



in the framework of AMDM models. 

We determine the parameters of the cosmological model which matches observa- 
tional data on the large scale structure of the Universe best and, marginalizing over 
all other parameters, we determine the upper limit for the neutrino content and the 
neutrino mass. We restrict ourselves the sub-class of models without tensor mode 
and neglect early reionization. 

EXPERIMENTAL DATA SET AND METHOD 

Our approach is based on the comparison of the observational data on the struc- 
ture of the Universe over a wide range of scales with theoretical predictions from the 
power spectrum of small (linear) density fluctuations. The form of the spectrum 
strongly depends on the cosmological parameters Q^, ^6, ^v, N^, h and Ug. If the 
amplitude on a given scale is fixed by some observational data then predictions for 
observations on all other scales can be calculated and compared with observational 
data. Minimizing quadratic differences between the theoretical and observational 
values divided by the observational errors, x^, determines the best-fit values of the 
above mentioned cosmological parameters. We use the following observational data 
set: 

1. The location £p = 197 ± 6 and amplitude = 69 ± 8fiK of the first acoustic 
peak in the angular power spectrum of the CMB temperature fluctuations deduced 
from CMB map obtained in the Boomerang experiment [2] (these data are sensitive 
to the amplitude and form of the initial power spectrum in the range ?a 200/i~^Mpc); 

2. The power spectrum of density fluctuations of Abell-ACO clusters 
{PA+Aco{kj)) obtained from their space distribution by [16] (scale range 10 — 
lOO/i-^Mpc); 

3. The constraint for the amplitude of the fluctuation power spectrum on ^ 
10/i~^Mpc scale derived from a recent optical determination of the mass function 
of nearby galaxy clusters [10]; 

4. The constraint for the amplitude of the fluctuation power spectrum on « 
10/i~^Mpc scale derived from the observed evolution of the galaxy cluster X-ray 
temperature distribution function between z = 0.05 and z — 0.32 [17]; 

5. The constraint for the amplitude of the fluctuation power spectrum on ^ 
10/i^^Mpc scale derived from the existence of three very massive clusters of galaxies 
observed so far at z > 0.5 [1]; 

6. The constraint on the amplitude of the linear power spectrum of density 
fluctuations in our vicinity obtained from the study of bulk flows of galaxies in 
sphere of radius 50/i^^Mpc V50 = (375 ± 85)km/s [13]; 

7. The constraint for the amplitude of the fluctuation power spectrum on 0.1 — 
lh~^Mpc scale and z — 3 derived from the Ly-a absorption lines seen in quasar 
spectra [11]; 

8. The constraints for the amplitude and inclination of the initial power spectrum 
on fti l/i~^Mpc scale and z = 2.5 scale derived from the Ly-a forest of quasar 



absorption lines [5]; 

9. Results from direct measurements of the Hubble constant h — 0.65 ± 0.10 
which is a compromise between by different groups; 

10. The nucleosynthesis constraint on the baryon density derived from a obser- 
vational content of inter galactic deuterium, fl^h'^ = 0.019 ± 0.0024 (95% C.L.) 
[4]. 

One of the main ingredients for the solution of our search problem is a reason- 
ably fast and accurate determination of the linear transfer function for dark matter 
clustering which depends on the cosmological parameters. We use accurate an- 
alytical approximations of the MDM transfer function T{k] z) depending on the 
parameters Qm, Q;,, and h given in Ref. [9]. The linear power spectrum 

of matter density fluctuations is P(/c; z) = A^ki^^T'^ik] z)Df{z)/Dl{0), where As is 
the normalization constant for scalar perturbations and Di{z) is the linear growth 
factor. We normalize the spectra to the 4-year COBE data which determine the 
amplitude of density perturbation at horizon scale, Sh [3]. Therefore, each model 
will match the COBE data by construction. The normalization constant ^4^ is then 
given by As = 27i^Sl{3000Mpc/hf+'''. 

The Abell-ACO power spectrum is related to the matter power spectrum at 
2; = by the cluster biasing parameter bd- PA+Aco{k) = bliP{k;0). We assume 
scale-independent linear bias as free parameter the best-fit value of which will be 
determined together with the other cosmological parameters. 

The dependence of the position and amplitude of the first acoustic peak in the 
CMB power spectrum on cosmological parameters Us, h, fib, flcdm and Oa can be 
determined using an analytical approximation given in [8] which has been extended 
to models with non-zero curvature (Jlfc = 1 — fl^ — Q,\ 0) in [7] . Its accuracy in 
the parameter range considered is better then 5%. 

The theoretical values of the other experimental constraints are calculated as 
described in [14]. There one can also find some tests of our method. 

RESULTS AND DISCUSSION 

The determination of the parameters 0,^, ^a, ^u, N^, ^6, h, Ug and hd by the 
Lcvcnbcrg-Marquardt minimization method is realized as follows: wc vary the 
set of parameters fi^, ^a, fly, fl^, h, Us and bd with fixed N,^ (= number of massive 
neutrino species) and find the minimum of x^. Since N^y is discrete we repeat this 
procedure for Ny—1, 2, and 3. The lowest of the three minima is the minimum of 
for the complete set of free parameters. We have seven continous free parameters. 
The formal number of data points is 24 but, as it was shown in [14], the 13 points 
of the cluster power spectrum can be described by just 3 degrees of freedom, so 
that the maximal number of truly independent measurements is 14. Therefore, the 
number of degrees of freedom for our search procedure is Np — -A^exp — A'par = 7. 
The model with one sort of massive neutrinos provides the best fit to the data, 
Xmm = Its parameters are following Vl^ = Q.STtoH, Qa = 0.69lo:2o, f^u = 



0.03l°:°L K = 1, = 0.037t°;lii, n, = 1.02^^0, h = O.Tll":?^, = 2At',l The 
errors are obtained by maximizing the (Gaussian) 68% confidence contours over all 
other parameters. 

However, there is only a marginal difference in xLin for iV^ = 1,2,3. With 
the given accuracy of the data we cannot conclude whether massive neutrinos are 
present at all, and if yes what number of degrees of freedom is favored. We summa- 
rize, that the considered observational data on LSS of the Universe can be explained 
by a AMDM inflationary model with approximately a scale invariant spectrum of 
scalar perturbations and small positive curvature, flk — —0.06. It is interesting 
to note that the values of the fundamental cosmo logical parameters Qm, and 
flk determined by these observations of large scale structure match the SNIa test 
Qrn — 0.75Qa = —0.25 ± 0.125 [15] very well. Models with vanishing A are outside 
of marginalized 3a contour of the best-fit model with = 1 determined by the 
LSS observational characteristics used here even without the SNIa constraint (Fig. 
la). For models with N^, = 2 and Ny = 3 contours are a bit wider, so models with 
vanishing A are outside the marginalized 2(t contour for arbitrary Ny. 

Results change only slightly if instead of the Boomerang data we use 
Boomerang-|- MAXIMA- 1. Hence, we can conclude that the LSS observational char- 
acteristics together with the Boomerang (-I-MAXIMA-1) data on the first acoustic 
peak already rule out zero-A models at more than 95% C.L. and actually demand a 
cosmological constant in the same bulk part as direct measurements. We consider 
this a non-trivial consistency check! 

The neutrino matter density Vly = 0.03 corresponds to a neutrino mass of niy ~ 
Mriyh'^ f» 1.4 eV but is compatible with within la. A ACDM model {Q^ = 0) is 
within the la contour of the best-fit AMDM model. 

To derive an upper limit for neutrino content and its rest mass, we first determine 
the marginalized la, 2a and 3a upper limits for fly for different values of Ny. Using 
the best-fit value for h at given Qy, we then obtain the a corresponding upper limit 
for the neutrino mass, niy = ^ATLyh? /Ny. The results are presented in Table 1 and 
Figure 1. For more species of massive neutrino the upper limit for Qy is somewhat 
higher but rUy is still lower for each C.L. The upper limit for fly raises with the 
confidence level as expected. But the upper limit for the mass grows only very 
little due to the reduction of the best-fit value for h. The upper limit for the 
combination flyh?' /N^'^^ is approximately constant for all numbers of species and 
confidence levels. The observational data set used here establishes an upper limit 
for the massive neutrino content of the universe which can be expressed in the form 
flyh^ /N^'^'^ < 0.042 at 2a confidence level. The corresponding upper limit on the 
neutrino mass rUy < 4eV is lower than the value obtained by [6]. 
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TABLE 1. Upper limits for the neutrino con- 
tent and mass (in eV) at different confidence 
levels. 





l<j C.L. 


2a C.L. 


3o- C.L. 








rrii, 


1 


0.10 3.65 


0.13 3.96 


0.18 4.04 


2 


0.15 2.79 


0.21 3.06 


0.29 3.35 


3 


0.20 2.40 


0.27 2.67 


0.35 2.78 




FIGURE 1. Likelihood contours (soHd Hne - 68.3%, dashed - 95.4%, dotted - 99.73%) of AMDM 
with = 1 in the fi^ ~ (a) and ~ (b) planes marginalized over all other parameters. 
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